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In a short superconducting nanowire connected to bulk superconducting leads, quantum phase 
slips behave as a system of linearly (as opposed to logarithmically) interacting charges. This system 
maps onto quantum mechanics of a particle in a periodic potential. We show that, while the state 
with a high density of phase slips is not a true insulator (a consequence of Josephson tunneling 
between the leads), for a range of parameters it behaves as such down to unobservably small tem- 
peratures. We also show that quantum phase slips give rise to multiple branches (bands) in the 
energy-current relation and to an interband ("exciton") mode. 



I. INTRODUCTION 

There is currently much interest in low-temperature 
properties of very thin superconducting wires. The key 
process, significance of which one aims to understand, 
is a quantum phase slip (QPS),— a virtual depletion of 
the superconducting density through which the system 
tunnels to a different value of the supercurrent. The 
rate of this process depends not only on the "bare" 
fugacity — the tunneling exponential associated with the 
phase slip core — but also on the interaction between in- 
dividual QPS. 

The present study has been motivated by the 
observation^ that, for a wire connected to bulk super- 
conducting leads, as the length of the wire is reduced, 
the logarithmic interaction between QPS crosses over 
to a linear one. Since we are talking about a tun- 
neling process, the interaction takes place in the [x, r) 
plane, where x is the spatial coordinate, and r is the Eu- 
clidean time. Heuristically, the change in the character of 
the interaction comes about because, during tunneling, 
the leads inhibit production of quasiparticles (plasmons) 
with nonzero wavenumbers, so the "lines of force" in the 
(x, t) plane align in the t direction and form a "string" 
with a nonzero tension. 

Now, it is well known that in the case of a logarithmic 
interaction, the system undergoes a phase transition that 
can be interpreted as unbinding of QPS-anti-QPS pairs. 
There are, in fact, several varieties of this phenomenon: 
a phase transition 3 - in the Kondo problem and the equiv- 
alent transition 4 ^ in the two-state dissipative quantum 
mechanics (DQM), the BKT transitional^ in the XY 
model and, finally, the superconductor-insulator transi- 
tion (SIT) in the DQM with a periodic potential.—^ Long 
superconducting wires with significant amounts of dis- 
order fall into the latter universality class ^ Is there a 
similar transition in short wires, where the interaction 
between QPS is linear? Experiments do bring out a dif- 
ference between shor t 12 ' 13 and lon g 14 ' 15 wires and suggest 
that a sharp transition exists in the former case. 

Here, we present two theoretical results concerning this 
problem. First, we point out that in a short wire connect- 
ing bulk superconductors there is no absolute distinction 
between a superconductor and an insulator: even an "in- 
sulator" can support a weak supercurrent due to Joseph- 



son tunneling. Second, we show that there is, neverthe- 
less, a range of parameters for which the wire acts as an 
insulator down to unobservably small temperatures. 

There is a mean-field-type approach 1 ^ to short wires 
that models QPS as an effective resistive environment. 
In that approach, the renormalization-group equation for 
the QPS fugacity is taken to be the same as in a long 
wire, except that the plasmon impedance is replaced by 
an effective shunt resistance. The latter is determined 
self-consistently and includes the QPS channel (R ps ), 
quasiparticles (R qp ), and the impedance of the electrodes 
(R e i ec ). While this procedure may qualitatively reflect 
the correct physics, it is not clear if it can be justified 
from first principles. Moreover, for the superconducting 
state in the limit when i? c ioc = and R qp 3> -R ps , it does 
not reproduce the exponential dependence of R ps on the 
temperature found 2 by direct calculation. 

Here, we consider the problem of SIT in short wires 
(with bulk superconducting leads) starting directly from 
the description in terms of a classical gas of particles 
with linear interactions in one (t) dimension (Sect. Hlf . 
We study this system in two different regimes (Sects. Hill 
and [Vj) , where two different approximations are possi- 
ble. In Sect. IIVI we pause to describe two potentially 
observable effects predicted by our theory: a breakdown 
voltage, which, under certain conditions, provides a di- 
rect measure of the QPS fugacity, and a curious "exciton" 
mode. Prospects for detection of the breakdown voltage 
in experiment are discussed in Sect. IVII Sect. I VIII is a 
conclusion. 



II. ONE-DIMENSIONAL GAS WITH LINEAR 
INTERACTIONS 

Our starting point is the partition sum of linearly in- 
teracting charges, which represent QPS and anti-QPS. 
This can be conveniently written as a path integral over 
an equivalent electrostatic potential </>(t): 

Z = J ^(r) e "*/ dr( ^ )2 x 

e i #wr- / n*^ 5 n ^e-^).(l) 

N±=0 +' ' J 2=1 m=l 
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Note that the integral over the constant (r-independent) 
component of <f> enforces the condition of charge neutral- 
ity: only terms with N + = N- contribute to the sum. 
Most of the expressions below will be for zero tempera- 
ture (i.e., infinite extent of the r dimension), but occa- 
sionally we will discuss consequences of the cutoff on r 
imposed by finite temperature. 

The parameters appearing in Eq. {1} are a, the bare 
QPS fugacity, and K, the strength of the linear interac- 
tion. The role of K can be understood by noting that, 
for a single QPS-anti-QPS pair, the path integral over 
4>{t) would be 

V(f)(T)e~^ I dT(d T 4>) 2 +i<p(Ti)-i4>(T2) _ g-|A"|ri-r 2 | (j> 

(2) 

These two parameters encode a wealth of microscopic 
details of the phase-slip process. For example, if QPS are 
rare, a is proportional to exp(— Score), where Score, the 
action of the QPS core, depends on the device specifics, 
such as the superconducting gap, the wire's width and 
thickness, and the width and thickness variations. When 
the latter are significant, QPS may favor a constriction, 
making a independent of the wire's length L. On the 
other hand, in a uniform wire, QPS occur preferentially 
at the center,— but the preference is rather weak, and a 
scales linearly with L. 

The value of K is determined by the one-dimensional 
superconducting stiffness K s {x) (including possible vari- 
ations of K s with the width and thickness): 



1 

K 



I 

4^2 



dx 
K 8 (x) 



(3) 



Here x is the coordinate along the wire. For a uniform 
wire, this gives K cx K s /L, i.e., K scales inversely with 
the length. The linear interaction in Eq. then matches 
that obtained in Ref. 0. We recall that, from the micro- 
scopic point of view, the linear interaction between QPS 
reflects the energy of the initial and final tunneling states, 
in which fluctuations of the magnitude of the order pa- 
rameter are small. This is why it is possible to find the 
strength of the interaction using the phase-only theory.^ 

In more detail, the relation of Eq. (fl]) to the phase-only 
theory is as follows. Consider the field <f>{x, r) dual to the 
phase 6{x,t) of the superconducting order parameter in 
the sense that d T 9 oc d x <j>. Since d x <j) is proportional to 
the charge density, <j> can be thought of as the density of 
the electric dipole moment (in suitable units). When the 
leads are bulk superconductors, all r-dependent compo- 
nents of 9 satisfy the Dirichlet boundary conditions (b.c.) 
at the ends of the wire^ In view of the duality relation 
above, the Dirichlet b.c. for 9 imply the Neumann b.c. 
for <f>, at both ends. 

Then, if a is much smaller than the charging energy 
of the wire, and the temperature is low enough, only the 
spatially uniform mode of 4>{x, r) matters, so instead of 
a field we have a single quantum-mechanical degree of 
freedom, 4>{t). Leaving aside overall constant factors, 



one can think of it as being the total dipole moment 
of the wire or, alternatively, the total charge transported 
through the wire since some fixed initial moment of time. 
It also coincides with the equivalent potential that ap- 
pears in Eq. JT]). 

Eq. {T]) is an effective low-energy theory, which encodes 
the high-energy details in the values of its parameters, 
and so is valid only up to some cutoff frequency A. Since 
the theory does not include fermionic quasiparticles, the 
cutoff is of the order of or smaller than the superconduct- 
ing gap A. Sometimes we will need a way to explicitly 
include such a cutoff in our calculations. One possibility, 
which we use in Sect. [Vj is to modify the kinetic term of 
so as to suppress modes above a cutoff frequency. An- 
other possibility, suitable for Monte Carlo simulations, 
is to discretize time. Cutoff dependence is discussed in 
Sect. IVI1 For now, we consider the case when both K 
and a are much smaller than A, so the cutoff does not 
matter. 

The sums in Eq. (fT|) can be evaluated by noting that 
the sum over N+ at fixed N = N + + AT_ is Newton's 
binomial, and the remaining sum over N is an exponen- 
tial (of a J dr cos <p). Thus, Eq. ([T]) is equivalent to the 
theory with the Euclidean action 



S 



(4) 



This is similar to the mapping between the Coulomb gas 
and the sine-Gordon model in two dimensionsi 17 i 18 ' 19 

Eq. (U)) is quantum mechanics of a particle in a cosine 
potential. Solutions to the corresponding Schrodingcr 
equation (the Mathieu functions) are known, and occa- 
sionally we will refer to their properties. For most of the 
time, however, we will consider various limits of Eq. ([J), 
for which the recourse to the exact solution is not neces- 
sary. 

Our system of charges with linear interactions is rem- 
iniscent of electrostatics in one dimension — an often 
used metaphor of quark confinement (see, for example, 
Ref. H3). As known in that context, depending on the 
dynamics of the charges, such a system can be in either 
a plasma or a confinement phase. The distinction lies in 
the large-distance behavior of a correlation function of 
fractional external charges. We therefore consider 



C,(r) = «e 



— 1 1 pil<t>f,T) -iq<)>{Si) 



» 



(•5) 



with arbitrary q. Double brackets denote averaging with 
exp(— S). In a plasma state, the charges can screen any 
external charge, integer or fractional. This corresponds 
to Eq. ([5]) approaching a constant value at large t, for 
any q. In the confinement phase, integer external charges 
can still be screened, by new charges nucleating from the 
vacuum, but fractional charges cannot. As a result, the 
correlator (JSJ) with fractional q goes to zero at large r. 

In the condensed-matter context, a plasma state, 
where the equivalent charges (the QPS) are unbound, 
corresponds to an insulator, and the confinement phase 
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to a superconductor. We now argue that, in the present 
case, there is no true insulator, i.e., Eq. ([5]) always goes 
to zero for fractional q. 



III. HIGH PHASE-SLIP DENSITY 

We begin with the limit that is the best candidate for 
an insulator: 



(6) 



when QPS are abundant. In this limit, the dominant 
paths are those where 4>(t) spends most of the time in 
the vicinity of a minimum of the cosine potential, but oc- 
casionally makes a transition to a neighboring minimum. 

A natural way to describe this physics is the semiclas- 
sical approximation. To take into account small fluctua- 
tions of <p near <fi = 0, we expand the cosine in Eq. Q in 
powers of 4> and treat the <f) 2 term as a "mass" , and the 
higher powers as an interaction. The Green function of 
</> in the free massive theory is 



where u> e is the characteristic frequency: 



ujt = Ka 



(7) 



(8) 



Using this Green function as a propagator for perturba- 
tion theory, one finds that the perturbative expansion is 
controlled by the small parameter K/a, and the asymp- 
totics of the correlator ([5|) at large r, in perturbation 
theory, remains close to 1. 

In addition to these perturbative fluctuations, however, 
there are large fluctuations (instantons) connecting dif- 
ferent minima of the cosine. They are solitons and anti- 
solitons of the sine-Gordon model, for example 



</>inst (t) = 4 arctan e u 



Their action is 



Si ns t = 8(a/A-) 1/2 



(9) 



These instantons are in a sense dual to the QPS, so when 
QPS are numerous, instantons are rare. Under the con- 
dition © , they form a dilute gas with average density 



inst ^ 



(10) 



Each instanton contained in the interval (0, r) con- 
tributes exp(27rig) to the correlator ([5]), and each anti- 
instanton exp(— lixiq). Statistics of these is that of an 
ideal gas, i.e., is given by the Poisson distribution P(M). 
Hence, 



C g (r) 



P(N+)P {Af- ) e 2 ™ q W+ - M - > 



-2rn.[l— cos(27rg)] 



We see that, due to the instantons, fractional charges 
are confined, with the string tension proportional to the 
instanton density. Note, however, that in practice r is 
cutoff by the inverse temperature, so the confinement 
of fractional charges (i.e., superconductivity) will not be 
observed until the temperature gets as low as T ~ n, 
which in the present case is exponentially small. 

The underlying physics becomes more transparent if 
one observes that dt4>/2w is the electric current in units 
of 2e (here t — —ir is the real time, and e is the electron 
charge). Thus, each instanton transports 2e of charge 
through the wire. A delocalized state of <f> (a Bloch wave 
of the theory @) corresponds to a steady current. At 
small K/a, the instanton rate is strongly suppressed, and 
for many purposes the system behaves as an insulator, 
but even then a weak steady current is possible by tun- 
neling of charges through the wire (the Josephson effect). 

Note also that, in a uniform wire, both K and a de- 
pend on the wire's length L, a being directly and K 
inversely proportional to it. Decreasing the length con- 
nects smoothly the "insulating" state of a nanowire to 
the superconducting state of a Josephson junction (JJ). 



IV. EXCITON AND THE BREAKDOWN 
VOLTAGE 

The spectrum of the theory Q has a band structure. 
In the path integral formalism, the quasimomcntum of 
these bands appears as a i) angle associated with the 
instanton number 



Q = 



l 

27 



drd T (, 



Thus, we generalize Eq. ((4]) as follows: 
S -» S + i$Q . 



(12) 



Each instanton now contributes an additional factor of 
, and in the dilute-gas approximation the dependence 
of the partition sum on $ is 



Z(0) 



exp[2/^(costf-l)] 



(13) 



The average current is obtained by differentiating In Z($) 
with respect to $ and equals 



I = tH^)) = 2nsini? 
2n 



(14) 



(11) 



where I is in units of 2e. The maximal (critical) current 
is I c = In. Eq. (fT4|) is of the familiar Josephson form, 
so we conclude that the quasimomentum (or •& angle) of 
the (^-description is precisely the relative phase A8 of the 
order parameter at the two ends of the wire. 

The interpretation of A9 as quasimomentum is partic- 
ularly apt given that in the absence of phase slips (i.e., 
for a — 0) A#, multiplied by the superfluid density, is 
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FIG. 1: A schematic of the band structure of the theory Q. 
The arrows show exciton transitions for two different values 
of the biasing current Note that the exciton frequency in 
general depends on I. 

the total momentum of the superfiuid in the wire. One 
can say that phase slips turn momentum into quasimo- 
mentum. This is similar to how phase slips enforce pe- 
riodicity with respect to the magnetic flux enclosed by a 
superconducting nanoringj^i 

Note, though, that lnZ(i)) gives only one band of the 
band structure, namely, the one corresponding to the 
ground state. In fact, there are additional bands, which 
can be obtained from the exact solution to the quantum- 
mechanical problem (H}. They correspond to additional 
branches of energy as a function of the quasimomcntum, 
as shown in Fig. [TJ One consequence of this is 
the existence of an interband excitation, which we term 
the exciton. At a fixed biasing current /, the exciton 
connects states for which the energy curves have equal 
slopes dE/dd = I, Since these states have different val- 
ues of d (in particular, for 1 = 0, the excited state is a 
7r-state), the transition produces a pulse of voltage. 

The exciton is distinct from the usual Josephson plas- 
mon, and the band structure described above is distinct 
from that described in Ref. I22I . The plasmon can be in- 
cluded in the theory by supplying i9 with nonstatic modes 
and generalizing the action (JT2J) into 

S = S+Y I drtidA + TY^P^K? , (15) 

where T is the temperature, f2 n = lixnT (n=integer) are 
the Matsubara frequencies, and P(Q) is the plasmon re- 
sponse kernel. The interaction between the plasmon and 
the exciton, represented by the second term in Eq. (|15p 
is the standard I'd interaction, with / cx d T 4> describing 
the fluctuating current. 

The only nonlincarity in Eq. (fT5|) is the cosine po- 
tential in S, Eq. ((4]). If i9(r) were the slowest variable 
in the problem, we could have, in the dilute instanton 
gas approximation, integrated (f> out, just as we did for 
a constant i? in Eq. (fT5)) . This is possible if, at small 



f2, P(£l) is purely capacitive, -P(fi) oc Cfi 2 , and the ca- 
pacitance C is sufficiently small. In this case, we re- 
cover the theory of Ref. [22- We know, however, that, in 
general, superconducting leads modify the low-frequency 
response significantly. For effectively one-dimensional 
leads, P(f2) oc \£l\/Z at small |f2|, Z being the plas- 
mon impedance of the leads (see Ref. d) , while for three- 
dimensional ones we expect -P(fi) to go to a constant. 
The limit of "bulk" leads considered in this paper corre- 
sponds to P being large, so that integrating out the n ^ 
modes of d results only in a small correction to the ki- 
netic term in Eq. This allows us to neglect plasmons 
altogether in a study of the low-frequency (ground-state) 
properties. We note, however, that the plasmon-exciton 
coupling may become important when the device is op- 
erated at a higher frequency, e.g., by being subjected to 
microwave radiation. 

Since <f> is proportional to the dipole moment of the 
wire, it couples, via a (j)V term, to the voltage V across 
the wire. For a wire deep in the "insulating" regime, cj> 
can sit for a long time near one minimum of the cosine 
potential. Its fluctuations there are small, and its value 
is related to the voltage by 

na 

V = — sin <p = V c sin <fi . (16) 

At small V , the exciton is a transition to an excited state 
near the same minimum. The critical value V c = ira/e 
is the breakdown voltage, at which <fr begins to escape 
classically, and the system becomes a conductor. 

Another breakdown channel, which opens at V > V g = 
2A/e, is production of quasiparticle pairs. However, the 
condition a <C A < A, under which we have been operat- 
ing (and validity of which for realistic wires is discussed 
in Sect. IVI[) . guarantees that V c <V g . 

V. LOW PHASE-SLIP DENSITY 

We now consider the limit of small QPS fugacity, 

a<min{if,A}, (17) 

where A is the ultraviolet (UV) cutoff. In this case, it is 
natural to expand the partition sum in powers of a. This 
brings in correlators of the Gaussian theory, 

So = ^kJ dT[d ^ )2 ' (18) 

which need to be regulated in the infrared. To this end, 
we modify Eq. (fl~8f as follows 

So^lJ dTdT'<f>(T)M (A, W r, r')^(r') , (19) 

where M$ is a differential operator whose inverse (the 
Green function) Go = M^ 1 is 

G (A, M ;r,r') = | (-e^ 1 ^' 1 \ e ~ hW ~ A ) • ( 20 ) 
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This implements an infrared (IR) cutoff at frequencies of 
order \x (,u <C A). Indeed, for \i\t - t'\ <C 1, Eq. (j2TJ|) 
nearly coincides with the unregulated Green function, 
proportional to r — r'| + const, but at large distances 
it decays exponentially, instead of growing linearly. In 
addition, Eq. (|20|) provides an explicit UV cutoff at fre- 
quencies of order A. 

The requisite correlators are 



e -i? 2 G (0) 



(21) 



(e ^(o) e T^M) = exp[-^ 2 + l)G (0)±gGo(r)],(22) 

etc. Single angular brackets denote averages in the Gaus- 
sian theory, and we have used the shorthand notation 
G (t) = G (A, li; t, 0). We do not include a $ angle in 
this section. 

For simplicity, in this section we calculate, instead of 
the correlator ((5]), directly its large r limit, i.e., the ex- 
pectation value of the "disorder parameter" 

((e < 9*(°)» = (e < «*W>+a / dr(e lqm cos c/>(t)) c + . . . , 



(23) 

where the subscript c denotes the connected part. The 
expectation value is even in q, and in what follows we 
assume q > 0. 

The first term on the right-hand side of Eq. (j2"3"| rapidly 
vanishes as the IR cutoff is removed (fx — * 0), but the 
second one generally does not. The relevant piece is 



dr(e 



iq<P(0) p -i<t>{r)\ 



qK 



'f(K/A,q) 



where / is a dimensionless function independent of fi; 
f(0,q) = 1. The approximation sign indicates that we 
have retained only the leading term in the li — > limit. 
We see that for any q ^ 1 Eq. (|2~4"1) vanishes exponentially 
at li — > 0, but for q = 1 it remains constant. 

The pattern repeats itself in higher orders in a, except 
that now larger (integer) charges can also be screened. 
This is again the physics of confinement. Note that the 
main contribution to the integral in Eq. (|24p for q = 1 
comes from small |t|, which means that the screening 
charge, represented by e~ l< ^ T ) in Eq. nucleates quite 
close to the external charge. For example, for K <C A, 
the main contribution comes from |t| ~ 1/K. This is 
how it should be since, according to Eq. the string 
tension is \K. 



VI. BREAKDOWN VOLTAGE IN 
EXPERIMENT 

At temperatures T <C A and with enough disorder, 
the ratio K/A depends on only one parameter — the total 
normal-state resistance of the wire Rn- 2 



K 
A 



2ir 2 R q 
Rn 



(25) 



Here R q = ir/2e 2 = 6.45 kfi is the resistance quantum. 
So, for a wire with R^ ~ R q , the ratio is K/A ~ 20. 
Estimates suggest that the right quantity for numerical 
comparison with the UV cutoff A < A is not K itself 
but jK or maybe even ^K. In either case, though, such 
a wire is quite far from the cutoff-independent regime 
K, a <C A. 

To locate the crossover between the "insulating" state 
and conventional superconductivity in the presence of a 
cutoff, we have used Monte Carlo simulations, with time 
discrctized in steps of At = tt/A and the crossover de- 
fined (somewhat arbitrarily) by the condition ((e 1 *)} = 
0.5. Preliminary data indicate that, at large K/A, the 
crossover occurs at a relatively small and nearly K- 
independent value a c w 0.4A. (Recall that a is the fugac- 
ity per the entire wire, and in a uniform wire is propor- 
tional to the length). Therefore, we do not exclude that 
the observed zero-bias peaks in differential resistances of 
"insulating" wire a 23 i 24 are a manifestation of the electri- 
cal breakdown effect described by Eq. (JTBJ) , and the ex- 
cess ( "offset" ) voltage contained in such a peak provides 
an estimate of the breakdown voltage V c . 



VII. CONCLUSION 

In this paper, we have presented a theory of short su- 
perconducting wires that uses, as its basic variable, the 
wire's electric dipole moment. We have found this de- 
scription to be well suited to cases when the more con- 
ventional phase variable has been made "heavy" (i.e., its 
fluctuations have become inhibited) by sufficiently large, 
"bulk" superconducting leads. The prominence of the 
dipole moment (polarization) means that we view a small 
superconductor as if it were almost an insulator. We have 
seen, however, that, even in a thinnest wire, the true insu- 
lating behavior is prevented by the Josephson tunneling 
between the leads. 

The thinner the wire is, the weaker that tunneling is. 
For sufficiently thin wires, it is described by instantons 
of the theory ^ , which are suppressed exponentially, so 
the insulating behavior persists down to exponentially 
low, unobservably small temperatures. The height of the 
potential traversed by the instantons, a in Eq. (f?|), has 
the interpretation of the bare fugacity of quantum phase 
slips in the wire. 

We have presented two potentially testable predictions 
of our theory. One effect, which, as discussed in Sect. IVI1 
may already have been observed, is the existence of a 
breakdown voltage in "insulating" wires. The other is the 
exciton, corresponding to the interband transitions in the 
theory It can presumably be detected as a feature 
in the absorption spectrum. We should note, however, 
that Eq. ^ can be literally applied to computation of 
the exciton spectrum only when both a and the exciton 
frequency are much smaller than the UV cutoff A. 

Note added in v4- Recently, we learned about the work 
of Mooij and Nazarov^ They start from a somewhat 
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different premise than we do, but their results overlap 
with ours. 
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